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ABSTRACT 

In a separate paper (Decin et al. 2003, henceforth paper I), we have re- 
examined the observations of IR excess obtained with the ISO satellite and dis- 
cussed the ages of stars with excess. The amount of dust (measured by the 
luminosity fraction fd = Lin/L+) seen around main-sequence stars of different 
ages shows several interesting trends. To discuss these results in the context of 
a physical model, we develop in this paper an analytical model for the dust pro- 
duction in Vega-type systems. Previously it has been claimed that a powerlaw 
slope of about -2 in the diagram plotting amount of dust versus time could be 
explained by a simple collisional cascade. We show that such a cascade in fact 
results in a powerlaw fd oc t -1 if the dust removal processes are dominated by 
collisions. A powerlaw fd oc t~ 2 only results when the dust removal processes 
become dominated by Pointing-Robertson drag. This may be the case in the 
Kuiper Belt of our own solar system, but it is certainly not the case in any of 
the observed disks. A steeper slope can, however, be created by including con- 
tinuous stirring into the models. We show that the existence of both young and 
old Vega-like systems with large amounts of dust (fd — 10~ 3 ) can be explained 
qualitatively by Kuiper-Belt-like structures with delayed stirring. Finally, the 
absence of young stars with intermediate amounts of dust may be due to the fact 
that stirring due to planet formation may not be active in young low-mass disks. 
The considerations in this paper support the picture of simultaneous stirring and 
dust production proposed by Kenyon & Bromley (2002b). 
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1. 



Introduction 



Debris disks are thought to form after the depletion of gas in young circumstellar disks. 
During the first ~ 10 Myr of a star's life, small dust grains grow by coagulation and finally 
produce planetesimals (Lissauer 1993). When the average velocity of the collisions increases 
due to stirring by planets or large planetesimals forming in the disk, the encounters between 
the planetesimals become destructive, initiating a collisional cascade. Large quantities of 
small dust grains are produced in such a cascade. These grains emit infrared and submm 
radiation, making the debris disk visible (e.g. Aumann et al. 1984; Zuckerman & Becklin 



Since most Vega-like excess stars are spatially unresolved by infrared telescopes, the 
quantity used to measure the amount of dust is the fractional luminosity, f d : the ratio of 
the dust emission Lir and the stellar luminosity 



fd can be interpreted as a covering fraction. It measures the fraction of the sky seen from 
the star which is covered by dust, and therefore the fraction of the stellar radiation which 
will be absorbed and reprocessed to the infrared. 

Studying the age dependence of the Vega phenomenon has been at the focus of several 
studies using the ISO satellite. The Vega phenomenon appears to be much more widespread 
around younger stars than around older stars. Studying a volume-limited sample of stars 
near the Sun, Habing et al. (1999, 2001) found that more than 60% of the observed stars 
below an age of 400 Myrs show a Vega-like excess, while only 9% of older stars can be 
classified as Vega-like. This result is also confirmed by submm observations (e.g. Zuckerman 
& Becklin 1993; Wyatt et al. 2003). Looking at a few selected vega-like stars, Kalas (2000) 
noted a ~ t^ 1 decline. Studying members of young clusters with different ages, Spangler 
et al. (2001) found a trend in the average amount of dust seen around stars of different ages 
which was described by a powerlaw dependence of fd versus time. Spangler et al. tentatively 
interpreted this result as a global trend due to a collisional cascade. 

A thorough study of the required parameters, age and fractional luminosity, was done 
in paper I to re-examine the available ISO data. It led to the following implications for 
models: (i) debris disks are more common around young stars than around old ones, but 
(ii) a general power-law for the dust mass versus age with slope of ~ —2, as found by e.g. 



1993). 
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Holland et al. (1998) and Spangler et al. (2001) could not be confirmed. Two main issues are 
causing the difference in results. There clearly is a strong decrease in circumstellar material 
from weak-line T Tauri stars to the youngest main sequence stars, but we reject weak-line 
T Tauri stars as debris disks (Lagrange et al. 2000) because these disks are likely to still be 
gas rich and governed by different physics (Artymowicz 1996). Also, we find a significant 
number of old stars with large fd- values, contradicting a global powerlaw decrease of fd- 
There may still be a decline of dust mass for the youngest stars, with slope of ~ -1.3, but 
surely not for the older ones, (in) We found that the maximum excess does not depend on 
age, also demonstrating the absence of a global declining trend applying to all stars. A final 
interesting result from paper I is that (iv) there is an apparent scarcity of young stars with 
intermediate or low infrared excesses. 

Little modeling on the general time evolution of Vega-like disks has been done so far. The 
collisional evolution of the solar system Asteroid Belt has been modeled extensively, mainly 
to understand the observed present day size and rotation distribution of large asteroids. 
Similar models for Kuiper Belt objects have been described by Davis & Farinella (e.g. 1997). 
Production of dust in the Asteroid Belt has been studied by Durda & Dermott (1997), 
focusing on the detailed size distribution produced by the collisional cascade at small sizes, 
which is influenced by both collisions and interaction with radiation. Studies of collisional 
cascades in the Kuiper Belt have also considered the formation of dust, but only in very basic 
ways (Stern & Colwell 1997b). The most detailed numerical models of combined growth and 
destruction of bodies in disks around stars are recent calculations by Kenyon & Luu (1999a,b) 
and Kenyon & Bromley (2001, 2002a,b) who study the evolution of the outer protoplanetary 
disk. We will get back to these papers in the discussion. 

The goal of the present paper is to better understand the age dependence of the dust 
mass in Vega-type stars in terms of a simple model. Therefore, an analytical model is devel- 
oped in Sect. 2 which describes the evolution of debris dust. This model is then confronted 
with the observations and discussed in Sect. 3. 



2. Collisional model 

We will use a very simple collisional model to derive the amount of dust produced by 
a collisional cascade of large bodies, and to compute the decay of the dust present in the 
system with time. A similar model has been used to estimate the lifetime of Vega-like disks 
by Kenyon & Bromley (2001), but we go further and compute the amount of dust as a 
function of time under the assumption that the evolution of the disk is purely due to a 
collisional cascade. In deriving the model, we will make use as much as possible of general 
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laws governing all kinds of collisions, and refer as little as possible to detailed material 
properties. The underlying model will be that of the Kuiper- or Asteroid belt of our own 
system: a number of large bodies, colliding to produce and replenish the dust we see. 



2.1. Collisional removal of comets 

Let us assume that all the dust in the system is produced as the end result of a collisional 
cascade involving large bodies. We therefore start with a number N c of comets. These comets 
have a radius a c , a geometrical cross- section for collisions between comets of <t co ii = 4na^, 
and a mass m c = ^-p c al, where p c is the density of the cometary material. 

For the collisions between the comets, we use a particle-in-a-box model: all comets are 
moving through a given volume V which could be an entire planetary disk or just a limited 
region like the Asteroid or Kuiper Belt. We assume that the occupied volume is a section of 
a wedge-shaped disk between distances r\ and r 2 with a local half-height H(r) = h-r where 
h is the normalized height H/r of the disk. The volume of such a wedge segment is 

V = y h{rl - r?) =: f hr* , (2) 

where f is a number close to r 2 if i\ is not too close to r 2 . If we use numbers typical for the 
Kuiper Belt (n ~ 35 AU, r 2 ~ 50 AU, h ~ .2) we find f = 43.5 AU and V = 6.25 x 10 44 cm 3 . 

The particle-in-a-box model assumes that there is a fixed collision velocity between 
different comets, which we denote with v co \\. For simplicity, we will write this collision 
velocity as a fraction v of the Kepler velocity: 



Veda = w k ( r ) = U V ~f~ ' ( 3 ) 

where G is the gravitational constant and M is the mass of the star. One can easily see 
that v is related to the relative height h of the disk. Collisions between bodies in a disk 
are due to differences in the inclinations and eccentricities of orbits. A particle on a circular 
orbit with inclination % will cross the midplane of the system with a vertical speed i>k x tani. 
Therefore, a disk with a normalized height of h = tani will show typical collision velocities 
h x f k . We can therefore assume v — h whenever we numerically evaluate expressions. 

With the collision velocity v co n and the volume given, we can compute the sweeping 
time t s which is the time needed by a single comet to sweep the entire volume accessible to 
the cometary cloud 

V 

t s = • (4) 

^collOcoll 
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Note that the sweeping time is actually independent of the collision velocity, since the volume 
occupied by the cometary cloud also grows with t> coll . The collision time for comets is then 
given by 

^cornet ~J^~ ■ (5) 

For dust production to occur, the collisions between comets must be destructive. As- 
suming that any such collision removes 2 comets from the cloud, we get for the time evolution 
of the number of comets 

dN c o iV c 2 



-2^ , (6) 



dt t 
with the solution 

N (t) = NM (7) 

For small times, the number of comets will be constant while for long times we find 

iV c (t)~- for t>^— . (8) 

The number of comets after a given time is therefore initially constant and later turns 
into a powerlaw with slope —1. 



2.2. Dust production 

Little is known about the dust production in collisions between comets or asteroids, as 
the dust production rate will depend to some extent on the internal structure of the body. 
If the body has already been fractured many times in previous collisions and now basically 
is a bag of sand, a massive collision will release large amounts of dust. If, on the other hand, 
the body still has some internal strength, it can be expected to produce a distribution of 
fragment sizes. While the experimentally determined distributions probable give good results 
for the larger fragments, the production rate of dust grains is very uncertain. Therefore 
we follow a different road and assume that dust production is the result of a collisional 
cascade where increasingly small particles collide with each other to produce smaller and 
smaller fragments. It has been known for many years that the size distribution produced 
by such a process assumes the shape of a powerlaw with the slope f{m) oc mT q or (for 
spherical particles) f(a) oc a 2 ~ 3q where q can be shown to be equal to 11/6 for a self-similar 
collisional cascade, largely independent of material properties (Dohnanyi 1968; Williams & 
Wetherill 1994; Tanaka et al. 1996). Significant deviations from this law appear only when 
the gravity component of the impact strength is considered - but for bodies smaller than 1km, 



- 6- 



the binding energy is dominating. The steady-state powerlaw is based on the competition 
between production and removal processes. At the smallest sizes, this equilibrium will be 
disturbed by the non-collisional grain removal processes like radiation pressure (blowout) 
and Poynting- Robertson drag (Durda & Dermott 1997). Also collisions with /5-meteorites 
may play a role (Artymowicz 1996; Krivova et al. 2000). The disturbance in the powerlaw 
can actually progress like a wave to larger particle sizes (Durda & Dermott 1997). For the 
simple estimates in this paper, we are going to ignore these additional processes. We will 
assume that a powerlaw distribution is produced by the collisions, and that the distribution 
continues with the same slope down to the particles which are removed from the system by 
another process. We write the size distribution as 

f(a) = W (9) 

where 7 = 2 — 3q is equal to —3.5. It is easy to show that the collisional timescales are 
slowest at the large-particle end of the distribution. Then, the collisions between large 
comets determines the injection of material in the cascade and therefore directly the dust 
production rate at the small particle end of the cascade. The steady state dust production 
rate TZ ga i n in such a model is proportional to the number of comet collisions per time, i.e. 

ft g ai„ = A;iiV c 2 , (10) 

where k± is a constant. 

Since the dust production rate is directly fixed by the collisional cascade, the dust 
removal processes govern the amount of dust visible in a given source. 



2.3. Dust removal 

Before we look into the details of dust removal, we can already provide the general time 
dependence of the dust content in a debris disk driven by a purely collisional evolution. From 
Eq. (10) it is clear that the dust production rate is proportional to rate of collisions between 
comets, thus oc iV c 2 . Since the number of comets decreases as i -1 , the dust production rate is 
proportional to t~ 2 . However, we cannot conclude from this that the amount of dust present 
in steady state shows the same time dependence - the mechanism which is responsible for 
grain losses has to be considered as well. If grain losses are due to collisions among the visible 
grains themselves, the loss rates will generally be proportional to the number of visible grains 
squared. However, if some other process removes the visible grains, the loss rates will only 
linearly depend on the number of the grains. Thus we have 



(11) 
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where a is either 1 or 2. In steady state we have 

j (1^) 

and therefore 

"-=S) 1/a ^ ■ < i3 » 

Thus, for internal collisional processes dominating the grain loss, we can expect the 
number of visible grains to be proportional to the number of comets left in the system, 
corresponding to a dependence of the dust amount. If other processes dominate grain 
removal, a iV c 2 (corresponding to t~ 2 ) dependence should be expected. 

In the following we will show under what circumstances which process dominates, and 
therefore which powerlaw we should expect from the observations. 



2.3.1. Poynting- Roberts on lifetime 

Poynting-Robertson drag removes angular momentum from particles in orbit around a 
star and causes them to spiral inward. The Poynting-Robertson drag reduces the size of a 
circular orbit according to (Burns et al. 1979) 

GM*2[3 2 (3 

r = — : — ~ = ~ 2u k- , (14) 



where c is the speed of light and (3 is the ratio of radiative and gravitational acceleration 
of a particle. The radiative force acting on the particle depends on the absorption and 
scattering cross sections (Burns et al. 1979). For the simple estimates in our study, we do 
not wish to discuss the detailed grain properties. Instead, we limit ourselves to large grains 
(compared to the wavelength of the radiation emitted by the star) and assume the grains 
to be perfect absorbers. In this limit, the absorption cross section of the particle is na 2 
and scattering does not contribute to the radiation pressure (van der Hulst 1981). (3 is then 
inversely proportional to the particle size, so we can write (Backman & Paresce 1993) 

0(a) = Z l± = 1^ (15) 

1 K ' l6ncGM*ap c 2 a v ; 

where is the blowout size, i.e. the size of particles which will be blown out of the system 
by radiation pressure. We define these to be the particles with (3 = 1/2. Particles with 
(3=1/2 will leave the system when ejected in a collision between particles unaffected by 
radiation pressure (Burns et al. 1979). 
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The PR timescale of particles of radius a moving in orbits with a semi-major-axis r is 
then given by 

T PR (a) = - = ^77^-— ( 16 ) 
r 2GM+ ab 

— V(^) 2 ^ ■ (17, 

Poynting-Robertson drag is frequently used in the literature to compute the dust loss 
in debris disks (e.g. Jura et al. 1998; Jourdain de Muizon et al. 1999; Li & Lunine 2003). It 
certainly provides an upper limit on removal time scales. However, as we will show below, 
collisions will always dominate the removal process in the currently observable debris disks. 



2.3.2. Collisional lifetime 

The collisional cross-section between two particles with radii ai and a 2 is given by 

a(a 1 , a 2 ) = n(a 1 + a 2 ) 2 . (18) 

Of course, not all collisions will be destructive - only collisions with particle sizes above a 
certain threshold can destroy the bigger particle. We denote by eai the size of the smallest 
particle which can (at a given collision velocity) still destroy a target particle with radius 
a±. The total destructive collision cross section provided by the particle size distribution in 
order to destroy particles with radius a\ is then given by 

0tot(ai) = J f(a)n(a 1 + afda (19) 

= J /a^Mai + afda . (20) 

In a powerlaw distribution with 7 < — 3, the cross section will be dominated by collisions 
with grain sizes eai. 

7re 7+1 

^tot(a) ^ -fa—ra^ =: e f a a^ (21) 
7 + 1 

where we have defined 

e = -7^+7(7 + 1) . (22) 
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The collisional lifetime due to collisions with other particles in the cascade is therefore 
given by 

V VI 
Tcoii(a) = rr = -r- t, (23) 

which for 7 = —3.5 is proportional to ^/a. However, the absolute value of the timescale 
depends on the normalization constant f a of the equilibrium size distribution. 

We still need to compute e. Collisions are destructive if the kinetic energy of the collision 
is approximately equal to the binding energy of the two bodies. In general, the binding energy 
is composed of a component due to the material properties, and a gravitational component. 
However, for bodies with 1 km or less in diameter, the gravitational component can be 
ignored. For a collision to be destructive we therefore require 

lvv 2 coll = S(m 1 + m 2 ) (24) 

where /x is the reduced mass and S is the binding energy per mass. We then ask what the 
minimum mass rri2 is which can still destroy a body with mass mi in a collision with a given 
velocity. Using 1712 = e 3 mi we find 



^3 ^coll . / V coU V Col\ 1 /ar\ 

6 " 4S ~Vl6^~^S ~ l ■ (25) 

Note that e contains the only important velocity dependence of our estimates. The higher 
the collision velocity, the larger the size interval of destructive impactors, and the larger the 
destruction rates of a given particle size. We will use a binding energy of S — 2 ■ 10 6 erg/s, 
a value typical for the icy bodies in the outer solar system (Kenyon & Luu 1999a). Table 1 
lists the value of eo as a function of S and i> co ii. 



2.3.3. Sublimation 

Sublimation of ice can be a removal process for small grains in the inner regions of the 
disk. Ice sublimation is an extremely strong function of the grain temperature. In a range 
between 90 and 120 K, the sublimation time r su bi can be approximated by the following 
expression (Backman & Paresce 1993) 

Tsubl ~ io 6 ^ ( -JjL_) ^ yr (26) 

where T g is the grain temperature. Vega-like disks show emission which is usually dominated 
by 60/im excess, only very few stars show detectable access at 25/im (e.g. Laureijs et al. 
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2002). The disk emission usually peaks around 60/im or at even longer wavelengths (e.g. 
Zuckerman & Becklin 1993). Using Wiens law, this corresponds to dust temperatures below 
100 K. The sublimation times for such temperatures are generally comparable with or longer 
than Poynting-Robertson drag timescales, even for the A stars among Vega-like stars (e.g. (3 
Pic, a Lyra and a Psa, see Backman & Paresce (1993)). In stars of later spectral type, grains 
at similar distance from the star will be colder. While grain sublimation may play a role in 
setting the inner boundaries of debris disks, it does not dominate the removal processes. We 
therefore ignore grain sublimation in the following. 



2.4. Steady-state size distribution 

In steady state in a collisional cascade, /„ can be determined equating the mass flux 
through the collisional cascade with influx of new material at the top end. The mass flux 
through the chain will in general depend on the details of the collision physics, but for the 
purpose of this study, a simple estimate is sufficient. The collisional mass loss of particles 
with sizes in a scale-free size interval between a and 2a is given by 

M(a) = /( " )a y = f a a ^ Wc£ ° . (27) 

Tcoll(o) V 

For the equilibrium slope 7 = —3.5, this value is independent of grain size. This mass loss is 
an excellent estimate for the mass flux through the cascade if in a typical collision the largest 
fragment has about half the diameter of the impactors. Equating M[a) with the inflow of 
mass into the collisional cascade (Eq. (6) multiplied by the m c ), we find 





f a = N c J 2 -^^ = N c a 5 J\l- . (28) 
-3-Pceo 



2.5. Dust visible in steady state 

The dust visible in the steady-state solution of a collisional cascade consists of two 
components. One component are small dust grains which are being blown away by stellar 
radiation (the wind grains). Each of these grains contributes to the radiation (scattering and 
IR/submm emission) only for the duration of about 1 Kepler time, typically a few hundred 
years. The second component consists of grains which are too large to be blown away, but are 
already small enough to provide significant surface for interaction with photons (the orbiting 
grains). Such grains contribute to the radiation as long as they live in the given orbit, i.e. 
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for a collision time or a PR-Drag time, whatever is shorter. The relative importance of both 
components can be easily estimated. The critical size for blowout of grains is a^. Of the grains 
smaller than this size, there will be a size a ra< \ which will dominate the radiation interaction 
of the grains below the blowout limit. We assume that the destruction of bigger particles 
produces the two different sizes according to the steady state powerlaw, i.e. the small 
grains are produced more frequently by factor (ab/a ra d)~ 7 - The production of grain surface 
(which is most important for the visibility of dust) proceeds then with a ratio (ab/a ra d)~ 7 ~ 2 . 
Therefore, the orbiting grains will dominate if their lifetime r meets the condition 

r > (a b /a rad ) _7_2 t ke p • (29) 

For a debris disk around a main-sequence star, the blowout size is typically 1 — 10/im (de- 
pending upon spectral type (Artymowicz 1988)) while a ra( j should be of the order of 1/im. 
Therefore, as long as the life time of the smallest orbiting grains is at least 10-100 orbits, 
orbiting grains will always dominate the visibility of debris disks. In the following discussion, 
we will therefore focus on the steady-state abundances of grains with sizes just above the 
blowout size. We will call this size of the visible dust grains a v ; s . 

2.6. Steady-state abundance of grains 

The number of grains visible in the disk at a given time will be determined by the 
dominating grain-loss mechanism. 

2.6.1. Collisionally- dominated grains 

If the visible dust grains are still being removed by collisions and have their steady-state 
abundance, we find 

n co \\ = f(a vis )da = f a a?da w N c a 5 J 2 J — a^ 1 (30) 

V eo 

where we have again assumed da ~ a. 

2. 6. 2. PR drag dominated grains 

If the visible grains are removed by PR drag, this means that the collisional cascade 
is effectively terminated at the size a vis . Then, the number of grains can be estimated by 
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equating the mass flux through the collisional cascade with PR driven losses. Therefore we 
have 

n PR = ^iV c rp R (a vis ) = iV c 2 ^^ . (31) 



2.6.3. General case 
The number of grains visible in the disk will then be given by 

n vis = min(n co ii,np R ) . (32) 
and the covering fraction /d of the disk is 

fd = n vis Q a bs(avis)7ra 2 is , (33) 

where Q a bs is the usual absorption efficiency of the dust particles. At high disk masses, the 
collisional destruction always dominates the life time of grains, while at low disk masses, 
PR drag can become important. The disk mass where the switch in processes takes place is 
given by the equation n co n = n PR . We find 

4tt /8tt y/a c a vis GM+r :i p c 

M C oii—PR = — \ / . (34) 

9 V eo c 



2.7. Dependencies and typical numbers 

The equations derived above show many properties of a collisional cascade and the 
amount of radiation reprocessing by dust which can be expected from such a cascade. If we 
insert typical numbers into the above equations, we find the following results: 
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^cornet ' • 4 

= 4.9 
Tcoll = 3.9 
tpr = 2.0 

Mcol^pR = 1.7 

/d,coii = 2.3 

where we have used a value of 226 for the e parameter which is appropriate for our standard 
model (S = 2 ■ 10 6 erg/g, v co n = 0.1v k ~ 0.45 km/s). 

First of all, we see that the dust removal timescales through collisions are much shorter 
than the removal timescales through PR drag. This means that typical debris disks will 
be collisionally dominated down to the smallest dust sizes which will then be blown out by 
radiation pressure. Eq. (39) shows that the disk mass (i.e. the mass of all comets in the 
collisional system) must be as low as 1.7- 10~ 3 M e for the transition to occur. From Eq. (40) 
we can see that the /d-value at that stage will typically be 4 • 10~ 7 - well below the detection 
limits of ISO and IRAS. 

Figure 1 shows the dependence of the dust visible in the system on the different pa- 
rameters of the model. Panel (a) shows the data taken from paper I. The solid line in the 
other five panels of the figure represents the time dependence of a standard model with 
Mdisk = 10M e , f = 43 AU, v = 0.1, a c = 1km around an AO star. The other lines indicate 
the changes due to the variation of a single parameter. From panel (b) in Fig. 1, it can 
be seen that only the more massive disks (> 10M e ) reach collisional steady state within 
Tcomet — 10 8 years. This is the time when most comets have seen at least one collision, and 
when the slope in the /d-time relation turns from zero to —1. Models with disk masses below 
10M e result in a constant dust production for 10 8 years or more, after which the curve turns 
into the common powerlaw with slope —1. 

The speed at which collisional equilibrium is achieved is also influenced by the size a c of the 
parent bodies (Fig le). If for a given disk mass, we reduce the size of the parent bodies, the 
number of such particles will increase and the collision time will become smaller. Collisional 
equilibrium is established after one or a few collision times. We can see this effect clearly in 
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the curves. For the standard size of 1 km, collisional equilibrium is reached after about 10 s 
years. For 10 km bodies this takes 10 9 years, and for 100 m bodies, collisional equilibrium is 
already fully established after 10 7 years. 

Another very important parameter is the collision velocity v = i> co ii/t>k (Fig. lc). Chang- 
ing the collision velocity changes the amount of dust seen in the system significantly, in a 
counter-intuitive way: increasing the collision velocity decreases the amount of dust seen 
from a steady-state cascade. This behavior results from the fact that the grain removal 
processes also strongly depend on the collisional velocities. Normally one would also expect 
that increasing the velocity would lead to faster depletion of the cometary cloud. However, 
this is not the case because the stirring of the orbits also increases the volume in which the 
particles move (see Eq. (2)) - the collision times of comets are therefore not influenced by 
the relative velocity. 

The distance of the cometary cloud (Fig lc) to the star also influences the observed values, 
in two ways. First, at a larger distance, the same amount of dust covers a smaller fraction 
of the solid angle seen from the star. The effect decreases the fd value proportionally to r~ 2 . 
Also, moving to larger distances increases all time scales. We can see in the figure that at 
larger distances, it takes much longer for collisional equilibrium to be established. The time 
of almost constant / d is extended to 10 9 years if we move the dust production site from 43 
to 150 AU. 

Finally, we can look at the dependence on the spectral type of the star (Fig. If). Moving 
from an AO star to later types slightly increases the amount of dust seen. The most impor- 
tant effect is that the size which dominates the visible dust will be smaller for low-luminosity 
stars since the blow-out limit moves to smaller grain sizes. An additional small effect is the 
lower mass of late-type stars which will reduce collision timescales. 

An additional parameter which influences the outcome of the calculations is the binding 
energy S, which is dependent on the internal composition and structure of the comets. The 
binding energy determines the value of e through Eqs. (22) and (25). This parameter influ- 
ences the speed at which material is processed through the collisional cascade and therefore 
the amount of visible dust. The calculation shown in this paper all use S — 2 • 10 6 erg/g, a 
value appropriate for icy comet-like bodies (Kenyon & Luu 1999a, and references therein), 
and we assume this value to hold for the entire size range in the collisional cascade. Using 
S = 10 7 erg/g, a value more appropriate for asteroid-type bodies, table 1 shows that the 
value of e would decrease by approximately a factor of 10. Eqs. (30), (36), (37), (39), and 
(40) all show a e^ 2 dependence, so the corresponding timescales and masses as well as 
the amount of visible dust would increase by a factor of about 3. Similarly, if the material 
the comets are made of would be exceptionally weak (Kenyon & Luu 1999a), the numerical 
values would decrease accordingly. 



- 15 - 



An important result is: nowhere in the parameter space covered by Fig. 1, a slope of 
t~ 2 is observed. This means that in all models shown here, the collisional removal of dust 
grains dominates over the PR-drag. This was to be expected from equation (39) where we 
showed that the transition from collisionally dominated dust removal to PR-drag dominated 
removal processes happens only at rather low disk masses, which are unobservable with 
current instrumentation. 



2.8. Summary of the collisional model 

1. Timescales do not depend critically on the collision velocities, if only the velocities are 
energetic enough to be destructive. 

2. A powerlaw dependence of the amount of observed dust as a function of time can 
only be expected after about one collisional time for the bodies starting the cascade. 
Before that time, an undisturbed cascade produces an approximately constant amount 
of dust. 

3. For lkm-sized comets, the collision times are of order of 10 9 years for disk masses 
(comets) of 1 Earth mass. To reach powerlaw behavior within 10 8 years and below, 
high disk masses (10M e ) are required. Alternatively, the collisional cascade could be 
started by 100 m bodies, provided that sufficient stirring can be achieved for these 
bodies. 

4. An undisturbed collisional cascade predicts /d = const for t < r comct and /d oc t^ 1 at 
later times. 

5. Disks in which PR drag dominates dust removal would show / d oc t" 2 behavior. How- 
ever, the transition from collisionally dominated to PR drag dominated disks happens 
at disk masses of typically 10 _3 M e , much less than required to support the observed 
debris disks. 

6. Adding more mass to the disk does not extend the lifetime of the debris disk. For 
given collisional velocities, disks of different masses all converge (after one collision 
time Tcomet) to the same curve. For more massive disks, this happens very quickly. 
The reason for this behavior is that in massive disks, the removal timescales are much 
shorter. In fact, in a more complete model which also treats planet formation, massive 
disks would form planets quickly and remove comets by gravitational scattering, further 
shortening the lifetime of the debris disk. 
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7. The lower blowout sizes of late type stars help somewhat to increase the amount of 
dust seen around old stars, but not enough to explain the old Vega-like stars with 
significant debris disks. 

3. Discussion 

In this section we will relate the findings of the simple model to the results of paper I, 
namely the main features of the distribution of stars in the log t - log fd diagram. 

3.1. The initial decrease for young stars 

Looking at only the youngest stars in the sample of paper I, and mainly at the A stars, 
there may indeed be an initial decrease in the amount of dust present. Due to the absence 
of very young stars with low /d values, it seems clear that the stars at the edge of the empty 
region in the lower left corner of the diagram must have evolved from larger fd values at 
younger ages. In this case, the lower boundary of the data points would be marking out the 
fastest path of decreasing the amount of dust in the system. Spangler et al. (2001) fitted a 
powerlaw with a slope —1.7 to their sample, but did include T Tauri stars which are probably 
not gas-free debris disks. For the reduced sample, the slope of the lower boundary is about 
— 1.3, much closer to the slope of —1 which we have derived for collisional cascades. The 
data clearly is currently not good enough to make a strong statement about the correctness 
of this slope. Hopefully, SIRTF will provide a much more solid database for this study. We 
would only like to make a remark here. Suppose the slope really is steeper than —1, what 
would that mean? A pure collisional cascade will not produce this. However, we can see 
from figure 1, that an increase in collision velocity is connected with a decrease in the fd 
value. Therefore, a slope steeper than —1 can be produced by a collisional cascade which is 
continuously stirred. We have simulated this in a very simple way, by increasing the collision 
velocity linearly in time starting from an initial v = 0.01 and ending at v = 1.0 after a given 
stirring time. The result of this experiment can be seen in Fig. 2. Compared to the curve 
with a fixed v — 0.1, the stirring curves all show a steeper slope. 

3.2. The upper limit 

Vega-like stars never seem to have fd values larger than 10 -3 . This upper limit is 
independent of age, i.e. at all ages do we find stars with fd ~ 10~ 3 , but not higher. This 
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upper limit can be understood from the cascade dynamics. Eq. (36) expresses the collision 
timescale for comets just like Eq. (35), but with the dependence on disk mass replaced by 
fd- From this equation we can see that the timescale for survival of the cometary cloud 
drops below 10 8 years when /d exceeds 10" 3 . An even stronger limit comes from work by 
Artymowicz (1996). He showed that in a gas-free debris disks, the amount of dust is limited 
by the creation of dust avalanches. When /d reaches values of 10~ 2 , the break-up of particles 
in the inner regions of the disk creates /9-meteorites which are small particles driven out by 
radiation pressure. For low f d values, such particles will not suffer additional collisions on 
the way out, because the disk is radially optically thin. However, at fd values between 1CT 3 
and 10~ 2 , the probability that such /3-Meteorites will cause additional collisions on their 
way out of the system becomes so high, that a self-accelerating avalanche effect is created. 
On basically a Kepler time, the disk will be cleaned of small visible dust grains. So even if 
the collisional equilibrium discussed in the present paper does allow larger fd values for a 
limited period of time, the avalanche mechanism effectively limits /d to about 10~ 3 . Stars 
with higher /d values must contain significant amounts of gas. The physics in such disks is 
different from "normal" Vega-like systems, and studying the age dependence must therefore 
make a clear distinction between the two types of disks. 



3.3. Independence of the upper limit of the age 

Models of collisional cascades indicate, that the dust content of a system should decrease 
with time. This seems to be inconsistent with the fact, that the observed debris disks can 
have / d ~ 1CT 3 at all ages. Assuming that the debris state starts at the same time for 
all stars, a systematic decrease of dust abundance should result in an absence of large /d 
values in old stars. This is not what is observed. We have seen above, that a more or less 
constant amount of dust can be supported before the disk goes into collisional equilibrium, 
i.e. if the parent bodies are large, if the distance from the star is large, or if the initial disk 
mass present in large planetesimals is small. However, figure 1 also shows that the types 
of solutions with constant dust for a long time always are at levels significantly lower than 
fd = 10 -3 . Therefore, a long duration of the debris state cannot be the explanation for 
old Vega-like stars. A more likely scenario is that different planetesimal disks are starting 
to become active debris disks at different times. We can study this in a toy model, in 
which we turn on the stirring of the planetesimal disk only after a given waiting period. 
Figure 3 shows the result of this calculation. In this way indeed solutions can be reached 
which cover the observed points. To find the reasons for the stirring, one needs to run 
much more complete models which self-consistently include the formation of planets. Such 
models have recently been introduced by Kenyon & Bromley (2001, 2002b). They run a code 
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which computes the time evolution of the size distribution in an annulus of a disk around a 
star. Processes included are growth through coagulation, collisional destruction, dynamical 
stirring and damping. It is found that stirring by bodies of 500 km or larger which form in the 
disk can lead to the onset of a collisional cascade which will produce debris dust (Kenyon 
& Bromley 2001). Furthermore, since the timescales for the growth of larger bodies are 
larger in the outer disk regions, planet formation and therefore debris production proceeds 
from the inner to the outer disk in rings and can lead to multiple collisional cascades at 
different times around the same star (Kenyon & Bromley 2002b). At least qualitatively, the 
ISO observations and the estimates put forward in the present paper are consistent with this 
picture. The question remains, if this mechanism can also explain very old Vega-like stars 
at ages of several Gyrs. 



3.4. The absence of young stars with low / d values 

The observations currently show no stars in this region. We have discussed in paper I, 
that the samples are currently to small to securely exclude the presence of stars in this 
region. However, lets for the moment take the relative absence of such stars as significant. 
Then this indicates than normal debris disk do not populate this region. Looking at figure 
1, there are three main ways to populate the area of young stars with low / d values. 

1. Low initial disk masses below 10M e . 

Our collisional model produces low / d values for low initial disk masses. However, the 
collisional velocity in the model is fixed. Considering a self-stirring scenario as proposed 
by Kenyon & Bromley (2002b), it is clear that enough mass must be present initially 
to allow for the formation of at least Pluto-like bodies. Kenyon & Bromley (2002b) 
use a model with 100M e of solid material to reach the required velocities within 100 
Myr. Stern & Colwell (1997a) find that at least 30M e are needed in the early Kuiper 
Belt in order to grow Pluto. Both calculations indicate that the time scale for planet 
formation and therefore for stirring is approximately inversely proportional to the disk 
mass. It is therefore likely that disks with less than about 10M e of solid material in 
the Kuiper Belt region will not be able to start a collisional cascade in the first 100 Myr 
because no large bodies will be formed in the belt. 

2. High collision velocities very early on 

Low /d values can also be produced by very high collision velocities in the disk, already 
very early, after about 10 7 years. However, such high velocities cannot be due to slow 
stirring by Pluto-sized objects forming in the disk. Large velocities are possible by 



embedding a giant planet in the disk, but in this case all planetesimals may have been 
removed by close encounters with the planet within the first 10 Myrs. 



3. Collisional cascade at large distance from star 

A collisional cascade at maybe 150 AU or further out seems to be capable of producing 
low fd values at young ages. However, in this case the question must be raised: if 
there is enough mass in the disk at such large radii to produce stirring and a collisional 
cascade, why is there not also material closer to the star which would form planets 
faster, start the collisional cascade earlier and produce dust more efficiently? If that 
additional material would be present, the star would become a "normal" debris disk 
with /d Rj 10 -3 . The collisional cascade at large distances will then only dominate the 
IR output of the system after the initial debris ring closer to the star looses its capacity 
of producing dust because most planetesimals in this regions have been removed. 

Therefore, we speculate that the absence of debris disks with intermediate or low / d values 
around young stars is due to the fact that stirring is only possible in higher mass disks. Low 
mass disks, if they exist, will remain quiet for the entire life of the star. 

4. Conclusions 

We have studied a simple model for collisional cascades in debris disks. While such a 
model does not lead to a self-consistent description of the stirring and dust production in 
a debris disks, it allows us to study the time behavior of such systems. Comparing with 
observations of debris disks, we come to the following conclusions. 

1. A collisional cascade with constant collision velocities leads to a powerlaw decrease 
of the amount of dust seen in a debris disk. The slope of that powerlaw is —1 for 
the parameter space valid for all observed debris disks. Only at much lower masses 
(typically 10~ 3 M e ) will this slope turn over to —2. 

2. A collisional cascade which is continuously stirred (i.e. where the collision velocities 
are increasing with time) can produce slopes steeper than — 1. 

3. If the initial decrease of fd for young stars is confirmed by further observations, it may 
be due to a combination of stirring and collisional cascade, as described by Kenyon & 
Bromley (2002b). 

4. The observed upper limit of fd ^ 10~ 3 for debris disks has to do with the dynamics of 
dust production in a collisional cascade and is due to the avalanche effect described by 
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Artymowicz (1996). According to that paper, larger /d values require gas to dominate 
the dust dynamics, but, following the definition by Lagrange et al. (2000), such disks 
are excluded from the class of debris disks. 

5. The most likely explanation for the presence of debris disk with / d values up to 10~ 3 
at ages well above a Gyr is the delayed onset of collisional cascades by late planet 
formation further away from the star. A prediction from this result is that the debris 
disks around older stars should be (on average) further away from the star than young 
debris disks. 

6. The tentatively observed absence of young debris disk with / d significantly lower than 
10~ 3 may be real, and caused by the effect of stirring. Low mass disks which could 
produce lower /d values cannot produce the planets needed to stir the disk quickly 
enough (Kenyon & Bromley 2002b). 

7. The trends observed in debris disks so far need confirmation by a much larger sample, 
hopefully available after the launch of SIRTF. 



We are grateful to Rens Waters for inspiring discussions on the subject. GD is supported 
by project IUAP P5/36 financed by the Belgian Federal Scientific Services (DWTC/SSTC). 
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Table 1: The e parameter as a function of collision velocity i> co u and binding energy S. 
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Fig. 1. — Dependence of the fractional luminosity of dust produced by a cloud of comets as 
a function of time. Panel (a): The observations (paper I). Explanation of the symbols: □: A 
main-sequence stars, A: F main-sequence stars, O: G dwarfs, x: K dwarfs and +: K giant. 
Panels (b)-(f) show the dependence on different parameters, starting from a standard model 
(solid line in all panels). See text. 
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Fig. 2. — /d-time relation with continuous stirring. Solid line: the standard model without 
stirring. The other lines include stirring, with v starting at 0.01 and increasing linearly with 
time to a maximum value of 1.0 after 10 8 , 10 9 , and 10 10 years. 
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Fig. 3. — /d-time relation for different starting times of stirring. Solid line: The standard 
model without stirring. The other lines show models with different starting times for stirring. 
The starting time (logt ) is noted at the curve. The stirring increases the velocity linearly 
between t and f • to where t is 10 7 , 10 8 , and 10 9 years. 



